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The physics of DPMJET-II.5 is described in a companion paper [Q. DPMJET samples hadron- 
hadron, hadron-nucleus, nucleus-nucleus and neutrino-nucleus interactions at high energies. The 
two-component Dual Parton Model is used with multiple soft chains and multiple minijets at each 
elementary interaction. Particle production is realized by the fragmentation of colorless parton- 
parton chains constructed from the quark content of the interacting hadrons. DPMJET-II.5 includes 
the cascading of secondaries within the target as well as projectile nuclei which is suppressed by 
the formation time concept. The excitation energy of the remaining target and projectile nuclei is 
calculated and using this nuclear evaporation is included into the model. At lab energies below 3 — 
5 GeV hadron-nucleus collisions are described within the conventional formation zone intranuclear 
cascade picture; thus the model may be applied down to such energies. It is possible to use the 
model up to primary energies of lO'^^ eV (per nucleon) in the lab. frame. 

DPMJET can also be applied to neutrino nucleus collisions. It extends the neutrino-nucleon 
models qel (quasi elastic neutrino interactions) and lepto (deep inelastic neutrino nucleon collisions) 
to neutrino collisions on nuclear targets. 
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I. PROGRAM SUMMARY 



Title of the program: 

Computer: 

Program language: 

Number of program lines: 

Other programs called: 
(included in DPMJET 
in modified form) 



Other programs called: 



DPMJET-II.5 

UNIX Workstations, LINUX-PC's 
FORTRAN-77 

about 80,000 (in addition the Lund codes linked) 
By DPMJET-II.5 

DIAGEN § 

Sampling of configurations for nucleus-nucleus interactions 

within the Glauber formalism 
DECAY |1 

Sampling the decay of hadron resonances. 
HADRIN g 

Sampling hadron-nucleon interactions below 5 GeV. 
parts of DTUJET [|§ 

Sampling of minijets and multiple soft chains, 
parts of FLUKA 

Nuclear evaporation and residual nucleus module. 
qel§ ^ ^ 

Quasi elastic neutrino-nucleon interaction. 
PYTHIA-6.1 § 

Sampling the hadronization of strings 

according to the Lund model |^ 

(Double precision version of JETSET ) 
lepto-6.5 |0) 

Deep inelastic neutrino-nucleon interaction. 

Note: lepto-6.5 uses the single precision JETSET-7.4. 



Method of solution: 



Monte Carlo event generator. 
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II. DESCRIPTION OF THE PROGRAM DPMJET-II.5 



The basic event generating unit of the code is the subroutine KKINC. Each call of this routine samples one inelastic 
hadron-nucleus or nucleus-nucleus interaction. Use and necessary initializations are described first in this section; a 
test program provided with the program package demonstrates the potential application. In the following subsection 
we discuss important model parameters and define their location in the code for potential user access. The basic 
structure of the supplied code is described in further subsections. 

A. The event generator and its initialization 

As already mentioned the code DPMJET-II.5 uses several other programs. The initialization of the DECAY |^ 
and HADRIN ||] codes requires one single call of the subroutines DPRIBL, DDATAR, DCHANT, DCHANH and 
DHADDE. 

Calls to PRBLM2 and JDTU initialize the multi-Pomeron sampling and the sampling of minijets like in the DTUJET- 
99 code (similar to DTUJET-93 (|§). 

A caU to PYINIT initiahzes the PYTHIA-6.1 samphng of chain decay. 

Further initializations and parameter definitions are provided in the routine DMINIT, via the named BLOCK DATA 
BLKDTl and the subroutine DEFAUL(EPN,PPN), which has to be called before event generation, too. Besides other 
initializations this routine sets the parameters, characterizing the type of the actual interaction, for vr+Cu collisions 
at 200 GeV; EPN and PPN (output variables) give energy and momentum in the lab-system. Actual predefinition 
and locations of those parameters which might be changed potentially by the user are discussed in the following 
subsection. 

Basic information for the application of the Glauber formalism according to the code DIAGEN [|j is generated by 
by the subroutine 

SHMAKI (IP, IPZ, IT, ITZ, RPROJ, RTARG, PPN) 
which requires the following input parameters defining the actual interaction: 

IP, IPZ: nucleon number/atomic number of the projectile nucleus; set IP=IPZ=1 for incident hadrons; 
IT, ITZ: nucleon number/atomic number of the target nucleus; 
PPN: projectile momentum in GeV/c (per nucleon) 

Since the calculations performed by SHMAKI are time consuming, in particular for heavy target nuclei, and in general 
have to be repeated for each different reaction type and energy, resp., the test program provided offers an option to 
prepare a data file 'GLAUBTAR.DAT' containing the necessary information for hadron-nucleus and nucleus-nucleus 
interactions to be considered. This file may be generated in a separate run of the test program using the option 
'GLAUBERF and GLAUBERA, compare Appendix A. The information from this file is read for a given projectile 
(projectile nucleus) (IP,IPZ) and target nucleus (IT, ITZ) by means of a subroutine call 

CALL SHMAKF(IP,IPZ,IT,ITZ). 
For the use of different target materials and or different projectiles in one calculation SHMAKF has to be called sub- 
sequently with the corresponding parameters (IP,IPZ,IT,ITZ). The information read from the file 'GLAUBTAR.DAT' 
is numbered internally by the index KKMAT=1,2,... according to the sequence of SHMAKF calls; up to 50 materials 
may be stored in the standard version of the program. If a larger number of materials is to be used the corresponding 
dimension in the common /DTUMAT/ has to be increased. (If the required data are not found in the file 'GLAUB- 
TAR.DAT' the execution is stopped.) 
After these initializations each call of the subroutine 

KKINC (EPN, IT, ITZ, IP, IPZ, IPROJ, KKMAT, ITARG, NHKKHl, IREJ) 
generates a single event. The input parameters not yet described have the following meaning: 

IPROJ: projectile type for hadron-nucleus collisions; 

DPMJET-II.5 uses the naming and internal numbering conventions from the BAMJET and 
DECAY [| codes which are listed in Table A— 2 of Appendix A.l; 

ITARG: Target type for hadron-hadron collisions; 

DPMJET-II.5 uses the naming and internal numbering conventions from the BAMJET |ll|,|l|] and 
DECAY [|| codes which are hsted in Table A-2 of Appendix A.l; 

KKMAT: controls the access of the event generator to the information on the Glauber formalism: 
KKMAT=0 : Glauber data expected from SHMAKI calculation; 

KKMAT>0 : Glauber data expected from the KKMAT'th caU of the subroutine SHMAKF. 
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NHKKHl: gives the position in the event COMMON HKKEVT, after which the final state particles 
are recorded. 



IRE J: IRE J = 1 indicates, that the event has been rejected 
IREJ = indicates, that the event is fine. 

If DPMJET-II.5 is used as event generator in a hadron cascade code, it is practical to write an interface routine. 
For the use of DPMJET-II.4 in the Cosmic Ray cascade code HEMAS-DPM (l|Jl|] there exists such an interface. 
DPMJET-II.4 has also been included into the CORSIKA code of Karhuhe 

The following commands will cause the generator to sample one inelastic tt+Cu event at 200 GeV laboratory energy: 

CALL DMINIT (NCASES , EPN , PPN , NCOUNT , IGLAUB) 
All of the following down to CALL SAMPPT is usually done within DMlNlT 

*** DECAY initialization 

CALL DDATAR 

CALL DCHANT 
*** HADRIN initialization 

CALL DHADDE 

CALL DCHANH 
*** setting default parameters 

CALL DEFAUL(EPN,PPN) 

CALL DEFAUX(EPN,PPN) 
*** JETSET initialization 

CALL LUND IN 

*** initialization of the random number generator supplied with DPMJET — II. 5 

CALL RNDMST(12,34,56,78) 
*** initialization for the Glauber formalism by explicit calculation 
*** for the actual reaction (projectile/target/energy = pi+/Cu/200 GeV) 

IP=1 

IPZ=1 

IT=64 

ITZ=29 

PPN=200 . 

CALL SHMAKI ( IP , IPZ , IT , ITZ , RPRD J , RTARG , PPN) 
*** sampling of 1 event (pi+ has type IPRDJ=13) 
IPRDJ=13 
KKMAT=0 

*** initialization of the evaporation module 

CALL BERTTP 

CALL INCINI 
*** initialization of the unitarization 

CALL PRBLM2(CMENER) 
*** initialization of the hard scattering 

CALL JTDTU(O) 

*** initialization of the transverse momenta for soft scattering 
CALL SAMPPT (0,PT) 



*** generating one event 

CALL KKINC (EPN , IT , ITZ , IP , IPZ , IPRD J , KKMAT) } 
STOP 

To use the information from the file 'GLAUBTAR.DAT' one has to call SHMAKF instead of SHMAKI and define 
KKMAT=1. 
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B. Information on the final state particles 



During the generation of single events several entries are scored in the common blocks /HKKEVT/ and /EXTEVT/ 
characterizing subsequent stages of the sampling process: Information on initial state nucleons as well as on interacting 
partons and constructed parton chains, decaying resonances and final state particles is stored into these common 
blocks. It has the following structure, completely defined in Appendix B: 

PARAMETER (NMXHKK= ) 

COMMON /HKKEVT/ NHKK.NEVHKK, ISTHKK(NMXHKK) , IDHKK(NMXHKK) , 
& JM0HKK(2,NMXHKK) , JDAHKK(2,NMXHKK) , 

& PHKK ( 5 , NMXHKK ) , VHKK ( 4 , NMXHKK ) , WHKK ( 4 , NMXHKK ) 

COMMON /EXTEVT/ IDRES (NMXHKK) , IDXRES (NMXHKK) ,NOBAM (NMXHKK) , 
& IDBAM (NMXHKK) ,IDCH (NMXHKK) ,NP0INT(10) 

The entries are characterized by their status ISTHKK and type IDHKK as well as by the 4-momenta PHKK; addi- 
tional pointers define 'parents' and 'daughters' of the actual entry. Final state particles are identified by their status 
lSTHKK(i) = l, -1 or 1001. Resonances, which have decayed are available with ISTIIKK(?;)=2. 



The structure of this common block closely follows the suggestions of Ref. [116 17 1; conventions for the description of 



the event history are described in some detail in Appendix B.2, followed by a sample event in Appendix B.3. 



C. Important parameters of the code 

In this subsection we summarize physical meaning and location of those parameters which may have a significant 
influence on observable quantities. 
Sampling of x values for partons 

Parton x-valucs for a given hadron are sampled from distributions of the general form 

q{x) (X x~°'{l — x)^ , 

subject to the reqirement that their sum is equal to one. The user has access to the /^-parameters for valence quarks 
which have, however, only minor influence on the final results. They are set to the following default values: 

(3^^'^ = 3.5, P^"" = 1.5; 

for reasons of an optimal sampling efficiency Psea = 0.0. From standard Regge arguments the power in x is fixed to 
be q;"'^" = 0.5 and a""' =0.5 for sea and valence quarks, respectively. 

Additionally, the x-values for partons of colliding projectile and target hadrons are correlated by the interaction 
mechanism within the DTU model: The constructed color-neutral parton-parton systems (chains) should aquire at 
least some minimum mass Mmm to allow the fragmentation into final state hadrons, 

^2 _ ^target . ^project. . ^ > j^2^_^. 

otherwise the collision numbers sampled according to the Glauber formalism may be strongly biased. Therefore, lower 
cuts are defined for the sampling of x-values, x™*" = C/-\/i, differing for valence quarks, diquarks and sea quarks, 
resp. To ensure minimum chain masses for valence-valence (v-v) systems the following values are set by default 
within DPMJET-11.5: 

C^"^' = 1.8, C™' = 2.0, C^"^^ = 0.5 

However, since these cuts are imposed for both the projectile and the target independently, lower (but still kine- 
matically allowed) chain masses tend to be suppressed at least for v-v systems. All the user accessible parameters 
discussed so far are located in the common block 

COMMON /XSEADI/ XSEACU,UNON,UNOM,UNDSEA, CVQ , CDQ , CSEA , SSMIMA , 
+ SSMIMQ.VVMTHR 

with the following assignments: 

= UNON P^"' = UNOM 

C™' = CVQ C™' = CDQ = CSEA 

M™™ = SSMIMA (Af™")2 = SSMIMQ M*^'' = VVMTHR 

Generation of transverse momenta for secondary hadrons 

In the model there are two sources of transverse momenta for created secondaries. First the partons aquire an 
internal p± which is taken into account in the construction of parton-parton chains. The subsequent hadronization 
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of the partonic chains is modelled by the JETSET code |9[] which itself assignes to the created particles transverse 
momenta with respect to the jet axis. 

Internal transverse momenta for partons are sampled within the subroutine SELPT (file dpmnucS.f) according to the 
following distribution for the reduced transverse energy Es = E±_ — nip, 

-^(xi^,exp(-72i^,/2) (2.1) 

with 

Pjl = y^Ef+2E~^, TTip = 0.94 GeV, 

where the slope parameter 7 = BBS is defined directly in this routine (default value BB3=6.0). However, in the 
case of severe kinematical limitations for individual partons/chains this parameter may be modified in the sampling 
process. 

Intranuclear cascade 

There are two important parameters controlling the development of the intranuclear cascade which have been 
defined in the previous section: The formation time tq of created secondary hadrons, and the factor a^^^ scaling the 
Fermi momenta of nucleons. With increasing tq the number of cascade generations and the number of low-energy 
particles will be reduced; a^^^ provides the possibility for a certain modification of the momentum distribution for 
low-energy nucleons. 

The maximum number of cascade generations KTAUGE (by default 0) has to be reset by the user, since the cascade 
is switched off by the default KTAUGE=0. 

In the code the corresponding parameters are defined in the common blocks 

COMMON /TAUFD/ TAUFDR , KTAUGE , ITAUVE 
COMMON /NUCIMP/ PRM0M(5,248) ,TAM0M(5,248) , . . . , 
+ PREBIN,TAEBIN,FERMDD,ETACOU 

where tq = TAUFOR and a^^^ = FERMOD. The option ITAUVE determines whether the pj_ -dependent definition 
of the formation time (ITAUVE=1) or the constant value tq (ITAUVE=2) are used; by default ITAUVE^l (compare 
Section 2.2). 

Lund JETSET fragmentation 

The chain fragmentation according to the Lund model JETSET contained within PYTIIIA-6.1 Q is customized 
for the DPMJET-II.5 needs by changing some of the Lund parameters in the initialization routine LUNDIN. More 
parameters are changed in the routine BAMLUN, where one JETSET fragmentation is called. 

In LUNDIN we prevent furthermore the weak decays of all otherwise stable hadrons and the decay of the tt" by setting 
the corresponding MDCY parameters equal to zero. 

Baryon and strangeness production, Popcorn mechanism: 

The popcorn mechanism in JETSET is controlled by the parameter pds (default: pdB ~ 0.1). 

COMMON /PDPCDR/ PDB.AJSDEF 

The parameters for the POPCORSE effect are in the COMMON /POPCCK/ and the parameters for the 
CASADIQU effect are in the COMMON /CASADI/. 

In JETSET pdB = 0. is used, to switch off the popcorn mechanism. 
Sea SU(3) symmetry: 

The parameter s'*'^" = SEASQ (default: SEASQ = 0.5) controls the s-quark content at the sea-quark chain ends. 

COMMON /SEASU3/ SEASQ 
Cronin efTect: 

The parameters MKCRO (default: MKCRO = 1) and CRONCO (default: CRONCO = 0.64) control the multiple 
scattering of partons within nuclear matter, and thus the Cronin effect. 

COMMON /CRONIN/ CRONCO , MKCRON 

MKCRON = switches off the Cronin effect. 

CRONCO is the parameter in the parton multiple scattering formula. 
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D. Structure of the supplied code 



The DPMJET-II.5 standard source code is stored in several FORTRAN files, listed in the following: 

• dpm25nucf .f, dpm25nuc2.f, dpm25nuc3.f, dpni25nuc4.f, dpni25nuc5.f, dpm25nuc6.f dpm25nuc7.f 

These modules contain all subroutines needed to describe the primary interaction according to the 
DPM as defined in the previous section. This includes the assignment of Fermi momenta to nucleons, 
construction of chains and monitoring their decay as well as the intranuclear cascade for generated 
secondaries. 

• dpm25hadri.f 

Includes the HADRIN Q routines for event generation in hadron-nuclcon collisions below 5 GeV as 
well as routines for the calculation of energy and reaction dependent cross sections. 

• dpm25nulib.f 

Contains few standard routines usually uneffected by further developments of the code, including the 
applied random number generator [ |l8| . 




• dpm25hist.f 

Contains a sample histogram routine producing line printer output for average multiplicities, rapidity 
and pseudorapidity distributions for generated particles. This serves mainly as an example on how to 
use the DPM JET events and for benchmarking the code when installing it. 

• dpm25diff.f 

Contains the routines used to call single diffractive particle production. 

• pythia61.f 

Contains double precision version of JETSET combined with the PYTHIA-6.1 code. 

• dpm25evap2.f or dpm25eva.f 

Contains the nuclear evaporation module. Please note, that the module dpm25evap2.f and the data 
file NUCLEAR.BIN can not be obtained from the author of DPMJET-II.5. The permission to use the 
evaporation module from FLUKA [Q has to be obtained directly from the Authors of FLUKA (Dr. 
A.Ferrari and Dr. P.R.Sala, CERN, Geneva and INFN, Sezione di Milano, 1-20133 Milano, Italy). The 
dpm25eva.f module contains dummy routines, to enable to run DPMJET-II without changes without 
the DPMEVAP module. Of course dpm25eva.f performs no evaporation of the residual nucleus. 

• dtu251ap.f or dpm251ap.f 

Contains the calculation of the minijets. 

Using the file dtu251ap.f the initialization of the minijets is done for one energy, given by the input 
cards. This version is used for stand-alone runs of the code, where each event is called for the same 
primary energy. 

Using the file dpm251ap.f this initialization is done for all energies. This version is suitable to insert 
DPMJET-II.5 into hadron cascade codes, where each event will be called at a new energy. 

• dtu25pom.f or dpm25pom.f 

Contains the calculation of multiple chain production according to eikonal unitarization. 

Using the file dtu25pom.f the initialization of the two component Dual Parton Model (calculation 

of the exclusive multi Pomeron distribution) is done for one energy, given by the input cards. This 

version is used for stand-alone runs of the code, where each event is called for the same primary 

energy. 

Using the file dpm25pom.f this initialization is done for all energies. This version is suitable to insert 
DPMJET-II.5 into hadron cascade codes, where each event will be called at a new energy. 
If the code is run repeatedly, it will be practical, to read these initialization tables, once calculated, 
from an external file. 
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• dpm25qelpo.f 

Contains the qel code ^ transformed to double precision and modified to use the Fermi momenta of 
the nucleons from DPMJET. 

• leptonew.f (and jetset74ku.f ) 

Contains the LEPTO-6.5 code in single precision, which uses the single precision jetset74.f. 

• dpm251epto.f 

Contains the DPMJET interface to LEPTO-6.5. 

• dpm25nonu.f 

If DPMJET is to be applied only to hadronic and nuclear collisions, not to neutrino collisions, it is 
possible to drop the files leptonew.f, jetset74ku.f, dpm251epto.f and dpm25qelpo.f and link the code 
instead to dpm25nonu.f, which contains only dummy routines. 

E. External input and output data files used 

• GLAUBTAR.DAT, Unit 47 

Glauber model data, see input cards GLAUBERI, GLAUBERA. 

• NUCLEAR.BIN 

Binary file to read in nuclear data for use in the evaporation module. 

• qel.evt 

Output file for the qel-events without the modifications of these events in DPMJET. 

• lepto.evt 

Output file for the lepto-events without the modifications of these events in DPMJET. 

F. Sequence of subroutine calls for event generation 

1. The generator subroutine KKINC monitors the sampling of a single event by one call of KKEVT, repeated 
FOZOKA calls and caU(s) of the subroutine FICONF. 

2. The subroutine KKEVT performs the sampling of the primary projectile-target interaction. This includes 

• sampling of the actual nucleon coordinates and assignment of 'partners' in the n elementary interactions 
between rip and rif nucleons from the projectile and the target, resp., done in SHMAKO 

• sampling of Fermi momenta of all nucleons in colliding nuclei in FER4M; 

• sampling of x-values for the appropriate quark systems (quarks, antiquarks and diquarks, resp.) from the 
interacting hadrons in the subroutine XKSAMP and assignment of quark flavors in FLKSAM; 

• construction of color neutral parton-parton chains in the subroutines KKEVVV, KKEVSV, KKEVVS, 
KKEVSS and further similar routines; 

• hadronization of chains via the JETSET code monitored by the subroutine HADRKK; 

• At energies below 3-5 GeV the code HADHAD is called as hadron-hadron event generator interfacing to 
the HADRIN code g. 

• Diffractive events are generated using a call to SDIFF. 

3. Alternatively the subroutines KKEVNU or KKEVLE treat the first step of the qel or lepto neutrino nucleus 
scattering. 

finally KKEVT, KKEVNU or KKEVLE returns the control to KKINC. 
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4. The intranuclear cascade is modeled by subsequent calls of FOZOKA for each secondary which did not leave 
the interacting nuclei. FOZOKA transports the particle either until the next interaction or until it leaves the 
nucleus. 

5. The Routine FICONF is called to perform the evaporation step. 



G. The test program DPMJET-II.5 



The event generator itself is supplied together with a test program which may serve as an example for the applica- 
tion of the program. By a call of the subroutine DMINIT (included in the source file dpm25nucl.f) the main program 
DPMJET-II.5 does the necessary initializations discussed above. It also allows the definition/modification of model 
parameters via input options explained in Appendix A. 

The sample histogram routine DISTR (source file dpm25hist.f) generates line printer output for average multiplic- 
ities, rapidity and pseudorapidity distributions for created particles. The scoring procedure demonstrates how the 
information on the event history stored in the common /HKKEVT/ may be applied to extract the properties of the 
final state particles. 



H. Remarks on the other codes applied in DPMJET— II.5: 
DIAGEN, DECAY, HADRIN and RNDM 



All the external codes applied are well documented. The generator DIAGEN |^ madified as described in samples 
configurations for nucleus-nucleus interactions in the framework of the Glauber model, i.e. spatial coordinates of 
nucleons in projectile and target nuclei, resp., as well as the actual numbers of 'elementary' interactions n, Up and rit 
(comp. sect. 2). DIAGEN routines are called via SHMAKI (initialization) and SHMAKO (sampling of configurations), 
both included in the file dpm25nuc7.f. 

A modified version of the resonance decay routine is defined by the subroutine DECHKK treating the decay of 
resonances from the common /HKKEVT/. 

The program HADRIN |^ for the description of hadron-nucleon interactions below 5 GeV is well tested against 
data [19). In DPMJET-II.5 it is called via the subroutine FHAD (file dpm25hadri.f). 

The algorithm applied for the generation of uniform random numbers (function RNDM) was described in ref. . 
The generated sequence of pseudorandom numbers should be independent of the actual hardware (if a real number 
has at least 24 significant bits in the internal representation). This generator has a period of about 2^'*'*. A few more 
comments and a test routine are provided in the source file dpm25nulib.f. 
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APPENDIX A : 
THE TEST PROGRAM DPMJET-II.5 



The test program demonstrates the application of the event generating routine KKINC and the extraction of 
information on the produced secondaries from the common block /HKKEVT/ (in the routine DISTR). Furthermore, 
it allows a simple redefinition of some important model parameters. It may also be used to prepare data files containing 
reaction specific information needed for the application of the Glauber formalism. In the case of event generation few 
standard histograms are constructed from sampled events. 

All program activities are monitored by input options. Each input option is identified by a code word and either 
changes the default values of variables and/or demands some action. Tab. 1 summarizes the standard input options. 
All variables read from the input file have default values. 



TABLE I. The most inijjortaiit oj)tions read by the snbroutiiic DMINIT 



Option 


Description 


TITLE 


next card is the run title 


COMMENT 


adds comments to the input data stream 


PROJPAR 


definition of the projectile 


TARPAR 


definition of the target 


GLAUBERI 


initialization of the Glauber formalism for 




hadron-nucleus interactions / data are written on unit 47 


GLAUBERA 


initialization of the Glauber formalism for 




nucleus-nucleus interactions / data are written on unit 47 


XSECNUC 


Calculation of Glauber cross sections 


HADRONIZ 


selects JETSET fragmentation of chains 


ENERGY 


beam energy {GeV; per nucleon for nuclei) 


MOMENTUM 


beam momentum {GeV/c; per nucleon for nuclei) 


CMENERGY 


cm. energy [GeV per nucleon) 


CENTRAL 


central A-A collisions forced 


TAUFOR 


definition of the formation time parameter 


SEADISTR 


monitors the x-behaviour of the quark distributions 


SEAQUARK 


monitors the x-behaviour of the sea-quark distributions 


SECINTER 


demands secondary interactions 


SEASU3 


defines s"*"" 


DIQUARKS 


selects sea-diquarks at chain ends 


POPCORN 


selects the popcorn mechanism 


POPCORCK 


selects the Capella-Kopeliovich popcorn mechanism 


POPCORSE 


selects the diquark breaking sea-quark mechanism 


CASADIQU 


selects Casado diagram 


CRONINPT 


selects the Cronin effect 


XCUTS 


monitors cuts for x-sampling 


FERMI 


monitors assignment of Fermi momenta to nucleons 


SINGDIFF 


controls the inclusion of single diffractivo events 


SINGLECH 


include Regge (single chain) contribution 


EVAPORAT 


include nuclear evaporation 


RANDOMIZ 


reinitialize the RNDM (RM48) randon number generator 


NOFINALE 


Skip routine FICONF, which calculates the residual nuclei 


START 


start event sampling for li-h, li-A or A-A collisions 


NEUTRINO 


start event sampling for qel neutrino collisions 


LEPTOEVT 


start event sampling for lepto neutrino collisions 


STOP 


stop the run 


PARTICLE 


print table of available particles and properties 


POMTABLE 


Write / Read a table for sampling multipomeron events 
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A.l: Description of input options and default parameters 



Table A 2 

Particle codes, SDUM is the parameter to be given on the PROJPAR input card. The internal codes and names 
were originally introduced by the BAMJET and DECAY [| codes. These codes are also used by the FLUKA 

§ and DTUJET |,§ codes.In the output of DPMJET, in COMMON HKKEVT, the particles are characterized with 
the PDG code, the internal code is given in the EXTEVT common. (Not all particles used internally are given in this 
Table.) 



oU U IVi 


internal code 


internal name 


particle 


PDG number 


charge 


ql 


^0 




r^xlU 1 UlN 


i 


r 


P 


zziz 


i 







z 


i 




o 
z 






-zziz 


-i 




-z 


o 
-z 


-i 




6 




e 


1 1 


-i 


n 
U 


n 
u 


n 
u 




4 




e 


1 1 
-ii 


i 


n 

u 


n 
U 


n 
U 




r 




ATT TT7 
IN U-Cj 




1 o 
iz 


U 


U 


U 


U 




D 


A ATT TTT 




-iz 


U 


U 


n 
U 


n 
U 




7 
/ 




7 


zz 


U 


U 


n 
U 


n 
U 


ATT7T TTT? AT 


Q 
O 


ATT7TT 


»,o 
n 


ziiz 


U 


o 
z 


i 


i 


A ATTTTTT'T^OAT 
AIN rj\J 1 xvUlN 


9 


A ATTTTT 


/Tin 


Olio 

-ziiz 


n 
U 


o 

-z 


-i 


-i 




iU 


IV/TTTTT 1 
iVl U rj+ 


,,+ 


-io 


1 
i 


A 

U 


n 
U 


n 
U 




ii 


1\/TTTT7 
M U tj- 




io 


-i 


U 


U 


U 




1 9 








n 


n 

u 



u 



u 


PION+ 


13 


PI+ 


7r+ 


211 


1 


2 


-1 





PION- 


14 


PI- 




-211 


-1 


1 


-2 





KAON+ 


15 


K+ 


K+ 


321 


1 


2 


-3 





KAON- 


16 


K- 


K- 


-321 


-1 


3 


-2 





LAMBDA 


17 


LAM 


A" 


3122 














ALAMBDA 


18 


ALAM 


AO 


-3122 














KAONSHRT 


19 


KOS 


Kl 


310 














SIGMA- 


20 


SIGM- 


n- 


3112 


-1 


1 


1 


3 
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SDUM 


internal code 


internal name 


particle 


PDG number 


charge 


ql 


q2 


q3 


SIGMA+ 


21 


SIGM+ 


S+ 


3222 


1 


2 


2 


3 


SIGMAZER 


22 


SIGMO 


SO 


3212 





2 


1 


3 


PIZERO 


23 


PIO 




111 





2 


-2 





KAONZERO 


24 


KO 




311 





1 


-3 





AKAONZER 


25 


AKO 




-311 





3 


-1 







31 


ETA550 


77 


221 





3 


-3 







32 


RHO+77 




213 


1 


2 


-1 







33 


RHO077 




113 





2 


-2 







34 


RHO-77 




-213 


-1 


1 


-2 







35 


OM0783 


LU 


223 
















36 


K*+892 


K*+ 


323 


1 


2 


-3 







37 


K*0892 


K*° 


313 





1 


-3 







38 


K*-892 


K*- 


-323 


-1 


3 


-2 







39 


AK*089 


K*° 


-313 





3 


-1 







40 


KA+125 


Kt 


10323 


1 













41 


KA0125 


K1 


10313 
















42 


KA-12r) 




-1032;-! 


-1 













43 


AKA012 


Ki° 


-10313 
















48 


S+1385 


S*+ 


3224 


1 













49 


S01385 


S*o 


3214 
















50 


S-1385 


E*- 


3114 


-1 













53 


N*++12 


A++ 


2224 


2 


2 


2 


2 




54 


N*+ 12 


A+ 


2214 


1 


2 


2 


1 




55 


N*012 


AO 


2114 





2 


1 


1 




56 


N*-12 


A^ 


1114 


-1 


1 


1 


1 




67 


AN~12 


A— 


-2224 


-2 


-2 


-2 


2 




68 


AN*-12 


A- 


-2214 


-1 


-2 


-2 


-1 




69 


AN*012 


AO 


-2114 





-2 


-1 


-1 
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SDUM 


internal code 


internal name 


particle 


PDG number 


charge 


ql 


q2 


q3 




70 


AN*+12 


A+ 


-1114 


1 


-1 


-1 


1 




95 


ETA* 




331 
















96 


PHI 


4> 


333 





3 


-3 







97 


TETAO 




3322 





2 


3 


3 




98 


TETA- 


2" 


3312 


-1 


1 


3 


3 




99 


ASIG- 




-3222 


-1 


-2 


-2 


-3 




100 


ASIGO 




-3212 





-2 


-1 


-3 




101 


ASIG+ 


E+ 


-3112 


1 


-1 


-1 


-3 




102 


ATETAO 




-3322 





-2 


-3 


3 




103 


ATETA+ 


S+ 


-3312 


1 


-1 


-3 


-3 




104 


SIG*+ 


E*+ 


3224 


1 


2 


2 


3 




105 


SIG*0 


S*° 


3214 





2 


1 


3 




106 


SIG*- 


s*- 


3114 


-1 


1 


1 


3 




107 


TETA*0 


77*0 


3324 





2 


3 


3 




108 


TETA* 


77* — 


3314 


-1 


1 


3 


3 




109 


OMEGA- 




3334 


-1 


3 


3 


3 




110 


ASIG*- 




-:-!224 


-1 


-2 


_2 


-3 




111 


ASIG*0 


E*° 


-3214 





-2 


-1 


-3 




112 


ASIG*+ 




-3114 


1 


-1 


-1 


-3 




113 


ATET*0 


— n 


-3324 





-2 


-3 


3 




114 


ATET*+ 




-3314 


1 


-1 


-3 


3 




115 


OMEGA+ 


n+ 


-3334 


1 


-3 


-3 


3 




116 


DO 




421 





4 


-2 







117 


D+ 


D+ 


411 


1 


4 


-1 







118 


D- 


D- 


-411 


-1 


1 


-4 







119 


ADO 


D" 


-421 





2 


-4 







120 


DS+ 


Dt 


431 


1 


4 


-3 







121 


DS- 


D7 


-431 


-1 


3 


-4 
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SDUM 


internal code 


internal name 


particle 


PDG number 


charge 


ql 


q2 


q3 




122 


ETAC 


Vc 


441 





4 


-4 







123 


D*0 




423 





4 


-2 







124 


D*+ 


D*+ 


413 


1 


4 


-1 







125 


D*- 


D*- 


-413 


-1 


1 


-4 







126 


AD*0 




-423 





2 


-4 







127 


DS*+ 


D*+ 


433 


1 


4 


-3 







128 


DS*- 


D*- 


-433 


-1 


3 


-4 







129 


CHIIC 


X\ 


20443 





4 


-4 







130 


JPSI 


J/^ 


443 
















131 


TAU+ 


T+ 


-15 


1 













132 


TAU- 


T~ 


15 


-1 













133 


NUET 


Vr 


16 
















134 


ANUET 


Vr 


-16 
















135 


NUEM 




14 
















136 


ANUEM 




-14 
















137 


LAMC+ 




4122 


1 


2 


1 


4 




i:-!8 


XIC+ 




42:-!2 


1 


2 


:-! 


4 




139 


XICO 


cO 


4132 





1 


3 


4 




140 


SIGC++ 




4222 


2 


2 


2 


4 




141 


SIGC+ 




4212 


1 













142 


SIGCO 




4112 





1 


1 


4 




143 


s+ 




4322 


1 













144 


so 


— /O 

'— 'c 


4312 
















145 


TO 




4332 





3 


3 


4 




146 


XU++ 


'— 'cc 


4422 


2 


2 


4 


4 




147 


XD+ 




4412 


1 


1 


4 


4 




148 


XS+ 




4432 


1 


3 


4 


4 




149 


ALAMC- 


K 


-4122 


1 


-2 


-1 


4 
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SDUM 


internal code 


internal name 


particle 


PDG number 


charge 


ql 


q2 


q3 




150 


AXIC- 


'c 


-4232 


1 


-2 


-3 


-4 




151 


AXICO 


'c 


-4132 





-1 


-3 


-4 




152 


ASIGC- 




-4222 


-2 


-2 


-2 


-4 




153 


ASIGC- 




-4212 


1 













154 


ASIGCO 




-4112 





-1 


-1 


-4 




155 


AS- 


'c 


-4322 


1 













156 


ASO 


'— 'c 


-4312 
















157 


ATO 




-4332 





-3 


-3 


-4 




158 


AXU- 


'—'cc 


-4422 


-2 


-2 


-4 


-4 




159 


AXD- 


^cc 


-4412 


1 


-1 


-4 


-4 




160 


AXS 


^cc 


-4432 


-1 


-3 


-4 


-4 




161 


Cl*++ 




4224 


2 


2 


2 


4 




162 


Cl*+ 




4214 


1 


2 


1 


4 




163 


C1*0 




4114 





1 


1 


4 




164 


s*+ 




4324 


1 


2 


3 


4 




165 


S*0 


77*0 

'c 


4314 





1 


3 


4 




1()() 


T*() 




4:-!:-!4 





3 


:-! 


4 




167 


xu*++ 


77*++ 
^cc 


4424 


2 


2 


4 


4 




168 


XD*+ 


77*+ 


4414 


1 


1 


4 


4 




169 


xs*+ 


Kt 


4434 


1 


3 


4 


4 




170 


TETA++ 




4444 


2 


4 


4 


4 




171 


ACl*- 




-4224 


2 


-2 


-2 


-4 




172 


ACl*- 




-4214 


-1 


-2 


-1 


-4 




173 


AC1*0 




-4114 





-1 


-1 


-4 




174 


AS*- 


'c 


-4324 


-1 


-2 


-3 


-4 




175 


AS*0 




-4314 





-1 


-3 


-4 




176 


AT*0 




-4334 





-3 


-3 


-4 




177 


AXU*- 


' — 'cc 


-4424 


-2 


-2 


-4 


-4 




178 


AXD*- 


^cc 


-4414 


1 


-1 


-4 


-4 




179 


AXS*- 


"cc 


-4434 


-1 


-3 


-4 


-4 




180 


ATET- 


^'ccc 


-4444 


-2 


-4 


-4 


-4 
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All input records of DPMJET-II.5 have the following form: 

CODEWD, (WHAT(I), I = 1,6), SDUM 
FORMAT (A8, 2X, 6E 10.0, A8) 

In the following we describe the meaning of the corresponding variables for the standard input options in the same 
order as the code words CODEWD are listed in Tab. A-1. 

Please note: further input options are described in the Appendices C and D. 

We give also the default values of the parameters . 

• Code word = 'TITLE' 

This option card must be followed by a card giving a run title, which will be reproduced in the output. 

• Code word = 'COMMENT' 

This option allows to add comments to the input file at arbitrary positions. 

WHAT(l): number of comment cards following this card, 
default : 1.0 

• Code word = 'PROJPAR' 

This card defines the type of the projectile; 

if given it has to be included before the MOMENTUM/ENERGY option(s). 

SDUM: defines the projectile to be a hadron if given; for naming conventions see Table 2. 

If SDUM is given WHAT(l) and WHAT(2) need no specification; for projectile nuclei SDUM has 
no meaning. 

WHAT(l): mass number of projectile nucleus - IP 
WHAT(2): atomic number of projectile nucleus - IPZ 
default: incident proton (IP=1, IPZ=1). 

• Code word = 'TARPAR' 

This card defines the type of the target nucleus. 
WHAT(l): mass number of projectile nucleus - IT 
WHAT(2): atomic number of projectile nucleus - ITZ 
default: Nitrogen N target ( IT=14, ITZ=7 ) 

• Code word = 'GLAUBERI' 

Requests the initialization of the Glauber formalism for hadron nucleus interactions; 
the target nucleus has to be defined by the code word TARPAR in advance; 

tables of impact parameter distributions for b-sampling are written to unit 47 for several momenta 

ipiab = •\/lO*'*^^''; « = 1, ...,24) and difi'erent projectiles (p, tt"*"); 

WHAT(l) = JGLAUB JGLAUB = 1 Calculation GLAUBTAR.DAT file (default) 

• Code word = 'GLAUBERA' 

Requests the initialization of the Glauber formalism for nucleus -nucleus interactions; 

the target nucleus has to be defined by the code word TARPAR in advance; 

the projectile nucleus has to be defined by the code word PROJPAR in advance; 

tables of impact parameter distributions for b-sampling are written to unit 47 for several momenta 

ipiab = a/Io''*'''^''; i — 1, 24) and different projectiles ; 

WHAT(l) = .JGLAUB JGLAUB = 1 Calculation GLAUBTAR.DAT file (default) 

• Code word = 'XSECNUC 
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Requests the calculation of Glaubercross sections for hadron (nucleus) nucleus interactions; 
the target nucleus has to be defined by the code word TARPAR in advance; 
the projectile nucleus has to be defined by the code word PROJPAR in advance; 
WHAT(l) = ECMUU lowest CMS energy 
WHAT(2) = ECMOO highest CMS energy 
WHAT(3) = NGRITT number of CMS energy points 

WHAT(4) = NEVTT number of Monte Carlo events for each cross section calculation 

Code word = 'HADRONIZ' 

Selects JETSET fragmentation of soft chains 

WHAT(l): IHADRZ (defauh: 2) 

IHADRZ = 2 selects JETSET fragmentation. 

IHADRZ = 11 selects an alternative JETSET fragmentation. 

Code word = 'ENERGY' 

This card defines the energy of the projectile in the target rest system. For incident nuclei the energy 
per nucleon is expected. 

NOTE: only one of the ENERGY and the MOMENTUM options is necessary, the last defined option 
is applied; both these options are to be given after the PROJPAR definition. 

WHAT(l): projectile energy in GeV 

Code word = 'MOMENTUM' 

This card defines the momentum of the projectile in the target rest system. For incident nuclei the 

momentum per nucleon is expected. 

NOTE: only one of the ENERGY and the MOMENTUM options is necessary, the last defined option 
is applied; both these options are to be given after the PROJPAR definition. 

WHAT(l): projectile momentum in GeV/c; 
Default: 100 000 GeV/c 

Code word = 'CMENERGY' 

Same as for code word 'ENERGY', but WHAT(l) defines the energy in the hadron/nucleon-nucleon 

cm. system. 

Code word = 'CENTRAL' 

This code word forces central nucleus-nucleus collisions, i.e. most nucleons of the projectile nucleus 
are forced to interact. The actual requirement depends on the atomic number of both the projectile 

and the target nuclei and is defined in the subroutine KKEVT (source file dpm25nuc2.f, after CALL 
SHMAKO ). Furthermore, the actual impact parameter is set near to zero for this case in subroutine 
MODE (source file dpm25nuc7.f). 

WHAT(l): = ICENTR default: 0.0, i.e. no forcing. 
ICENTR=0 : normal collisions 

ICENTR=1 : central collisions with impact parameter condition in dpm25nuc7.f and NA condi- 
tion in dpm25nuc2.f 

ICENTR=2 : only NA condition in dpm25nuc2.f 
ICENTR=3 : less central condition for Pb-Pb 
ICENTR=10: peripheral collisions 

Please note: The definition of central collisions is not unique, see the actual definitions in routine 
KKEVT (dpm25nuc2.f and dpm25nuc7.f). 

Code word = 'TAUFOR' 

This option defines the formation time parameter controlling the intranuclear cascade. Additionally 
it allows to restrict the number of generations of secondary cascade interactions. 
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WHAT(l): formation time in fm/c\ 
default: 105 /m/c 

WHAT(2): maximum number of allowed generations of secondary interactions; default: 
WHAT(3): monitors the definition of the formation time actually applied (comp. Subsects. 2.3 and 
3.3.3): 

WHAT(3)=1 : p_L -dependent formation time , (default) 
WHAT(3)=2 : constant formation time 

Please note, for h-h interactions, one should use WHAT(2) = 0. ! The default corresponds to a 
supression of the formation zone cascade. To apply the formation zone cascade, the recommended 
values are (TAUFOR 5. 25. 1.). 

Code word = 'SEADISTR' 

This option card defines properties of the quark distributions, which are of the general form q{x) oc 

WHAT(l): no meaning, a^^" is now controlled by the SEAQUARK card, (default a''"' = 0.5) 

default : 0!^*=" = 1.0 
WHAT(2): /3""'^ for valence-quark distributions of nucleons; 

default : = 3.5 
WHAT(3): (3^^^ for valence-quark distributions of mesons; 

default : =1.11 
WHAT (4): no meaning; 

NOTE: for reasons of the sampling efficiency the parameters for the sea distribution arc fixed to 
a*'^"=either 1.0 or 0.5 and p^^"' = 0.0 in the present version of the subroutine XKSAMP and cannot 
be changed easily by the user. 

Code word = 'SEAQUARK' 

This option card defines properties of the quark distributions, which are of the general form q{x) oc 

WHAT(l): a^^"" for sea-quark distributions possible values: 0.5 and 1.; 

default : a'"'' = 0.5 

NOTE: for reasons of the sampling efficiency the parameters for the sea distribution are fixed to 
a*'^"=either 1.0 or 0.5 and /3**" = 0.0 in the present version of the subroutine XKSAMP and cannot 
be changed easily by the user. 

Code word = 'SECINTER' 

This option demands secondary interactions, the use of this option is at present only recommended 
for heavy ion collisions at SPS energies. 

WHAT(l): = 1. secondary interactions demanded 

default: 

SECINTER should only be used for heavy ion collisions at CERN-SPS energies. 
Code word = 'SEASU3' 

This option determines the strange quark fraction s*®" at the sea-quark chain ends. 

WHAT(l): 8^*=", 
default: 0.5 

Code word = 'DIQUARKS' 

This option determines diquarks at sea chain ends 

WHAT(l): IDIQUA, 

default: 1 IDIQUA=0 no diquarks at Glauber sea chain ends 
IDIQUA=1 diquarks at Glauber sea chain ends 
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WHAT(2): IDIQUU, 

default: 1 IDIQUU=0 no diquarks at unitary sea chain ends 
IDIQUU=1 diquarks at unitary sea chain ends 

WHAT(3): AMEDD, 

default: 0.9 1.- AMEDD is the fraction of sea chain ends with a sea-diquark 

Code word = 'POPCORN' 

This option determines the popcorn mechanism for JETSET fragmentation. 

WHAT(l): PDB, 
default: 0.10 

JETSET: PDB gives the fraction of diquarks fragmenting directly into baryons. 

PDB = switches off the popcorn mechanism. 

NOTE: POPCORN should appear before the HADRONIZE card. 
Code word = 'POPCORCK' 

This option determines the CK (Capella-Kopeliovich) popcorn mechanism for JETSET fragmentation. 

WHAT(l): IJPOCK, 
default: 

IJPOCK=0 switches off the CK popcorn mechanism 
mechanism. 
WHAT(2): PDBCK, 

default: 0.00 

PDBCK gives the fraction of diquarks with CK mechanism 

NOTE: POPCORCK should appear before the HADRONIZE card. 
Code word = 'POPCORSE' 

This option determines the CK (Capella-Kopeliovich) popcorn mechanism for JETSET fragmentation. 

WHAT(2): PDBSE, 
default: 0.45 

PDBSE gives the probability for a diquark to be split by Glauber sea quarks 

WHAT(2): PDBSEU, 
default: 0.45 

PDBSEU gives the probability for a diquark to be split by unitary sea quarks 

NOTE: POPCORSE should appear before the HADRONIZE card. 
Code word 'CASADIQU' 

This option determines the CASADO diagram 

WHAT(l): ICASAD, 
default: 1 

ICASAD=0 switches off the Casado mechanism 
ICASAD=1 switches on the Casado mechanism 

WHAT(2): CASAXX, 
default: 0.50 

CASAXX gives the probability for the Casado mechanism 
Code word = 'CRONINPT' 
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This option determines the Cronin effect. 

WHAT(l): MKCRON, 
default: 1. 

MKCRON = switches off the Cronin effect. 
WHAT(2): CRONCO, 
default: 0.64 

CRONCO is the parameter in the parton multiple scattering formula. 
Code word = 'XCUTS' 

This option redefines the lower cuts for the sampling of x-values to ensure minimum chain masses for 
hadronization (used in XKSAMP). 

WHAT(l) = CVQ ; (x™')™" = CVQ/V^, 

lower cut for valence quarks; default : CVQ = 1.8 

WHAT(2) = CDQ ; (x™')™" = CDQ/^i, 

lower cut for valence diquarks; default : CDQ = 2.0 
WHAT(3) = CSEA; (a;^^p™" = CSEA/^/i, 

lower cut for sea quarks; default : CSEA = 0.5 

WHAT(4) = SSMIMA; (^xseayarget . (^sepproject. . ^ > {SSMIMAf, 

lower cut for the mass of sea-sea chains applied in XKSAMP; 
defauh : SSMIMA = 1.2 GeV. 

Code word = 'FERMF 

WHAT(l): Inclusion of Fermi momenta for nucleons if WHAT(1)= 1.0 (default) 

WHAT(2): FERMOD - scale factor for Fermi momentum as calculated from zero temperature Fermi 

distribution of nucleons; 

default: FERMOD = 0.6. 
WHAT(3): use zero temperature Fermi momentum distribution if WHAT(3)= 1.0 (default); use 

distribution according to Ref. [|o) for WHAT(3)=2. 

Code word = 'SINGDIFF' 

This option controls the generation of single diffractive events. 



WHAT(l): 


ISINGD; 


ISINGD 


= 1: 


Single diffraction included, 


ISINGD 


=0: 


Single diffraction supressed, 


ISINGD 


=2: 


Only single diffraction, 


ISINGD 


=3: 


Only single diffraction, target excited, 


ISINGD 


=4: 


Only single diffraction, projectile excited 


default: 


1 





Code word = 'SINGLECH' 

WHAT(l): ISICHA; 

ISINGD=1: include Regge (single chain) contributions, 
ISINGD=0: single chains supressed, 
ISINGD=2: Only single contribution, 
default: 

Please note: The single chain (Regge) contributions are only essential at low energies, they are at 
present only implemented for antibaryon and meson projectiles. 

Code word = 'EVAPORAT' 

Evaporation is performed if the EVAPORAT card is present Default:IEVAP=0 No evaporation 

what (1) = lEVAP 

lEVAP — 1: evaporation is performed 
lEVAP = 0: no evaporation 
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Actually the EVAPORAT card allows for more options, which normally are not needed by the user. 
See thee description in dpm25nucl.f. 

Code word = 'RANDOMIZ' 

Re-initialization of RNDM (RM48) random number generator Default: standard initialization 
ISEEDl and ISEED2 are printed at the end of each DPMJET run. They can also be generated by a 
RD20UT call (AUAUAU=RD20UT(ISEED1,ISEED2)) 

what (1) = ISEEDl 
what (2) = ISEED2 

Code word = 'NOFINALE' 

This option skips the call of the routine FICONF, for instance in heavy ion collisions; 
WHAT(I) = 1. : no FICONF call, Default: 0.: 

Code word = 'START' 

This option starts the generation of events in h-h, h-A and A-A collisions including the output of 

standard histograms. 

WHAT(l): number of events to be sampled; default: 100 
WHAT(2): the Glauber initialization is forced to be calculated in SHMAKI, 
i.e. no data read from file GLAUBTAR.DAT, if WHAT(2) = 1.0 

Code word = 'NEUTRINO' 

This option starts the generation of events for quasi elastic (qel) neutrino-nucleus collisions including 
the output of standard histograms. 

WHAT(l): number of events to be sampled; default: 100 

WHAT(2): neutrino type as defined in qel (1=1^6, 2=Pe, 3=1^^, 4=P^, 5=i'r, &=Dr) 

Code word = 'LEPTOEVT' 

This option starts the generation of events for deep inelastic (lepto) neutrino-nucleus collisions in- 
cluding the output of standard histograms. 

WHAT(l): number of events to be sampled; default: 100 

WHAT(2): neutrino type as defined in lepto {12=Ue, -12=9^-, 14=1'^, -14=i/^) 

Code word = 'STOP' 

This option stops the execution of the program. 

Code word = 'PARTICLE' 

This card triggers a printout of all the particles defined in the BAMJET-DECAY chain fragmentation, 
including name conventions, quantum numbers and decay channels. This are the internal particle 
codes used in DPMJET. 

Code word = TOAITABLE^ 

Only if the files dpm25pom.f and dpm251ap.f are linked instead of the files dtu25pom.f and dtu251ap.f. 
WHAT(l): ( default) the file pomtab.dat is written. :1 The file pomtab.dat is read. 
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A. 2: Sample Input 



In the following several typical examples of input data are reproduced. We give only the input cards which change 

the default parameters. 

One of the examples demonstrates the use of the test program to generate the input data for the Glauber formalism 
(file GLAUBTAR.DAT) in A-N colUsions. 

1. Hadron-Hadron collisions 
Note: SINGDIFF includes diffractive events. XCUTS slightly changes some defaults. 

TITLE 

DPMJET p-p 200 GeV LAB w. Diffr. 
PROJPAR PROTON 
TARPAR 1.0 1 . 

MOMENTUM 200 . 
SINGDIFF 1 . 

XCUTS 0.70 2.0 0.30 1.201 

START 1000000 . 1 . 

STOP 



Hadron-Nucleus collisions 



Note: Here the NOFINALE and EVAPORATE options assure evaporation and the working out of the residual 

nucleus. If this is required, then also the TAUFOR card is needed in a form as given to perform the FZIC. In 
hadron-nucleus collisions diffractive events might be demanded like in the example or might be supressed. 



TITLE 

DPMJET p-N 100000 GeV CMS w. Diffr. 
PROJPAR PROTON 
TARPAR 14.0 7. 

CMENERGY 100000. 
NOFINALE 0. 
EVAPORATE 1 . 

TAUFOR 5.0 25 . 

SINGDIFF 1 . 

START 50000 . 1 . 

STOP 



Minimum bias nucleus-nucleus collisions 



Note: Here again the NOFINALE and EVAPORATE options assure evaporation and the working out of the residual 
nuclei. If this is required, then also the TAUFOR card is needed in a form as given to perform the FZIC. In nucleus- 
nucleus collisions diffractive events might be demanded or might be supressed like in the example. In very high energy 
nucleus-nucleus collisions the code might run more stable if the Cronin effect is suppressed like in the excample given. 



TITLE 

DPMJET Fe-N 10 TeV CMS without Diffr. 

PROJPAR 56. 26. 

TARPAR 14.0 7. 
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CMEMERGY 10000. 

NOFINALE . 

EVAPORATE 1 . 

TAUFOR 5 . 25 . 

SINGDIFF . 

CRONINPT 0. . 00 

START 10000 . 1 . 

STOP 



Central Nucleus-Nucleus collisions 



Note: In central micleus-micleus collisions difFractivc events should always be excluded like in the examples given. 
In both examples the FZIC and evaporation is demanded (TAUFOR, NOFINALE, EVAPORATE options) It is also 
possible to run nucleus-nucleus collisions without FZIC and evaporation. The last parameter on the XCUTS card is 
convenient to tune somewhat the secondary particle multiplicity in nucleus-nucleus collisions. 



TITLE 

DPMJET central S-S 200GeV LAB with secondary interactions 



PROJPAR PROTON 

PROJPAR 32.0 16. 

TARPAR 32 . 16 . 

MOMENTUM 200 . 

NOFINALE 0. 

EVAPORATE 1 . 

TAUFOR 5.0 25 . 

SEC INTER 1 . 
CENTRAL 2 . 

SINGDIFF . 

CRONINPT 0. 0.00 

XCUTS 1.80 2.0 0.50 2.50 

START 10000. 1. 

STOP 



TITLE 

DPMJET central Pb-Pb 158 GeV LAB with secondary interaction 



PROJPAR 


207.0 


82. 


TARPAR 


207.0 


82. 


MOMENTUM 


158. 




NOFINALE 


0. 




EVAPORATE 


1. 




TAUFOR 


5.0 


25. 


SECINTER 


1. 




CENTRAL 


2. 




SINGDIFF 


0. 




CRONINPT 


1. 


0.64 


XCUTS 


1.00 


2.0 


START 


1000. 


0. 


STOP 
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Gross section calculation 



TITLE 

DPMJET p-Ar XSEC 
PROJPAR 

TARPAR 40 . 

XSECNUC 2 . 

STOP 



PROTON 



18. 
1000000. 



50. 



2000. 



TITLE 

DPMJET Fe-N XSEC 
PROJPAR 



PROJPAR 
TARPAR 
XSECNUC 
STOP 



56.0 
14.0 
2. 



PROTON 



26. 
7. 

1000000. 



50. 



20000. 



Input data for Glauber formalism 



TITLE 

DPMJET-II Prepare Glauber data for different nucleus-N collisions 



PROJPAR 1 . 1 . 

TARPAR 14.0 7. 
GLAUBERI 

PROJPAR 4. 2. 

TARPAR 14.0 7. 
GLAUBERA 

PROJPAR 12. 6. 

TARPAR 14.0 7. 
GLAUBERA 

PROJPAR 14. 7. 

TARPAR 14.0 7 . 
GLAUBERA 

PROJPAR 24. 12. 

TARPAR 14.0 7. 
GLAUBERA 

PROJPAR 28. 14. 

TARPAR 14.0 7 . 
GLAUBERA 

PROJPAR 40 . 20 . 

TARPAR 14.0 7. 

GLAUBERA 

PROJPAR 48 . 22 . 

TARPAR 14.0 7. 

GLAUBERA 

PROJPAR 56.0 26. 

TARPAR 14.0 7. 
GLAUBERA 
STOP 
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qel neutrino interactions in DPMJET-II.5 



Please note: DPMJET does not consider neutrinos as projectile particles, therefore, we use PRO JPAR = PROTON. 
This is of course not used later on. 



TITLE 

DPMJET-II.5 5 GeV Neutrino (nu-e)-Ar quasi-elastic interaction 

PROTON 

18. 



25. 
1. 



PRO JPAR 




TARPAR 


40.0 


MOMENTUM 


5. 


NOFINALE 


0. 


EVAPORATE 


1. 


TAUFOR 


5.0 


NEUTRINO 


1000. 


STOP 





Please note: DPMJET does not consider neutrinos as projectile particles, therefore, we use PRO JPAR = PROTON. 
This is of course not used later on. 



Deep inelastic (lepto) neutrino-nucleus interaction in DPMJET-IL5 



TITLE 






DPMJET Deep 


inelastic 


Neutrino (nu- 


PRO JPAR 






TARPAR 


40.0 


18. 


MOMENTUM 


50. 




NOFINALE 


0. 




EVAPORATE 


1. 




TAUFOR 


5.0 


25. 


LEPTOEVT 


1000. 


12. 


STOP 







PROTON 
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APPENDIX B : 

EVENT HISTORY AND THE COMMON BLOCKS /HKKEVT/ AND /EXTEVT/ 

B.l: Structure of the common blocks 

During the generation of individual events several entries are scored in the common blocks /HKKEVT/ and /EX- 
TEVT/ characterizing subsequent stages of the sampling process. Scored entries are, for instance, initial state nucle- 
ons, partons and parton chains, decaying resonances as well as final state particles. These entries are characterized 
by their type, 4-momenta and coordinates; additional pointers define 'parents' and 'daughters' of the actual entry 
(if any) . The structure of the /HKKEVT / common block closely follows the suggestions of Ref . ||l6|,|l^ . Within the 
code there are extensive comments explaining the variables used in this common block. Below the common blocks 
are reproduced together with these comments. 

Note that for interactions potentially resulting in high multiplicities of secondaries (i.e. very high energies and/or 
heavy ion-ion collisions) it may become necessary to increase the dimension NMXHKK, the maximum number of 
entries for a given event. 

PARAMETER (NMXHKK=49998) 

COMMON /HKKEVT/ NHKK , NEVHKK , ISTHKK (NMXHKK) , IDHKK (NMXHKK) , 
& JMOHKK (2, NMXHKK) .JDAHKK (2, NMXHKK) , 

& PHKK ( 5 , NMXHKK ) , VHKK ( 4 , NMXHKK ) , WHKK ( 4 , NMXHKK ) 

COMMON /EXTEVT/ IDRES (NMXHKK) , IDXRES (NMXHKK) , NOBAM (NMXHKK) , 
& IDBAM (NMXHKK) ,IDCH (NMXHKK) ,NPDINT(10) 

C 

C Based on the proposed standard COMMON block (Sjostrand Memo 17.3,89) 
C 

C NMXHKK: maximum numbers of entries (partons/particles) that can be 

C stored in the common block. 

C 

C NHKK: the actual number of entries stored in current event. These are 
C found in the first NHKK positions of the respective arrays below. 
C Index IHKK, 1 <= IHKK <= NHKK, is used below to denote a given 
C entry . 
C 

C ISTHKK ( IHKK) : status code for entry IHKK, with following mesinings: 



C = : null entry . 

C = 1 : an existing entry, which has not decayed or fragmented. 
C This is the main class of entries which represents the 

C "final state" given by the generator. 

C = 2 : an entry which has decayed or fragmented and therefore 
C is not appearing in the final state, but is retained for 

C event history information. 

C = 3 : a documentation line, defined separately from the event 
C history, (incoming reacting 

C particles, etc.) 

C = 4 - 10 : undefined, but reserved for future standards. 

C = 11 - 20 : at the disposal of each model builder for constructs 
C specific to his program, but equivalent to a null line in the 

C context of any other program. One example is the cone defining 

C vector of HERWIG, another cluster or event axes of the JETSET 

C analysis routines. 

C = 21 - : at the disposal of users, in particular for event tracking 

C in the detector. 

C 



C IDHKK(IHKK) : particle identity, according to the Particle Data Group 

C standard. 

C 

C JMOHKK (1 , IHKK) : pointer to the position where the mother is stored. 
C The value is for initial entries. 
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c 

C JM0HKK(2,IHKK) : pointer to position of second mother. Normally only 
C one mother exist, in which case the value is used. In cluster 
C fragmentation models, the two mothers would correspond to the q 
C and qbar which join to form a cluster. In string fragmentation, 
C the two mothers of a particle produced in the fragmentation would 
C be the two endpoints of the string (with the range in between 
C implied) . 
C 

C JDAHKK(1 , IHKK) : pointer to the position of the first daughter. If an 
C entry has not decayed, this is 0. 

C 

C JDAHKK(2 , IHKK) : pointer to the position of the last daughter. If Ein 
C entry has not decayed, this is 0. It is assumed that the daughters 
C of a particle (or cluster or string) are stored sequentially, so 
C that the whole range JDAHKK(1 , IHKK) - JDAHKK(2, IHKK) contains 
C daughters. Even in cases where only one daughter is defined (e.g. 
C KO -> KOS) both values should be defined, to make for a uniform 
C approach in terms of loop constructions. 
C 

C PHKK(1,IHKK) : momentum in the x direction, in GeV/c. 
C PHKK(2,IHKK) : momentum in the y direction, in GeV/c. 
C PHKK(3,IHKK) : momentum in the z direction, in GeV/c. 
C PHKK(4,IHKK) : energy, in GeV. 

C PHKK(5,IHKK) : mass, in GeV/c**2. For spacelike partons, it is allowed 
C to use a negative mass, according to PHKK(5,IHKK) = -sqrt (-m**2) . 
C 

C VHKK(1,IHKK) : production vertex x position, in mm. 

C VHKK(2,IHKK) : production vertex y position, in mm. 

C VHKK(3,IHKK) : production vertex z position, in mm. 

C VHKK(4,IHKK) : production time, in mm/c (= 3.33*10**(-12) s) . 

C 

C WHKK(I,IHKK) gives positions and times in projectile frame, 

C the chains are created on the positions of the projectile nucleons 

C in the projectile frame (target nucleons in target frame); 

C both positions are threfore not completely consistent. 

C The times in the projectile frame , however, are obtained by 

C a Lorentz transf ormtion from the lab system. 

C 

C Entries to the /EXTEVT/ COMMON block 

C 

C IDRES(IHKK) gives the Mass number A in case of nuclear 
C fragments or residual nuclei (IDHKK( IHKK) =80000) 
C 

C IDXRES(IHKK) gives the charge Z in case of nuclear 
C fragments or residual nuclei (IDHKK( IHKK) =80000) 
C 

C NOBAM(IHKK) = 1 for particles resulting from projectile fragmentation 
C =2 for particles resulting from target fragmentation 

C 

C IDBAM(IHKK) gives the internal dpmjet particle code (BAMJET code) 
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B.2: Conventions for the scoring of the event history in /HKKEVT/ 

In the following we briefly characterize the subsequent entries to /HKKEVT/ together with the most important 
conventions for their classification. 

• projectile hadron/nucleons: 

for projectile nuclei Fermi momenta in the projectile rest frame and coordinates within the nucleus are stored 
in the arrays PHKK and VHKK, resp., by the subroutine KKEVT; 

interacting and non-interacting nucleons have the status ISTHKK=11 and 13, resp.; 
nucleons wounded by the formation zone intranuclear cascade get status 17. 

• target nucleons; 

Fermi momenta in the target rest frame and coordinates within the nucleus are defined in PHKK and VHKK, 

resp. (KKEVT): 

interacting and non-interacting nucleons have the status ISTHKK=12 and 14, resp.; 
nucleons wounded by the formation zone intranuclear cascade get status 18. 

• valence quarks / diquarks from the interacting projectile hadron/nucleon(s) defined in the subroutine XKSAMP 

(total number IXPV); 

PHKK(3,...)=PHKK(4,...) contains the actual momentum fraction, VHKK the position of the 'mother' hadron; 
defined status ISTHKK=21; 

• sea quarks from interacting projectile hadrons defined in XKSAMP (total number IXPS); 
PHKK(3,...)=PHKK(4,...) contains the actual momentum fraction, VHKK the position of the 'mother' hadron; 

defined status ISTHKK=31; 

• valence quarks / diquarks from interacting target nucleons defined in XKSAMP (number IXTV); 
PHKK(3,...)=PHKK(4,...) contains the actual momentum fraction, VHKK the position of the 'mother' hadron; 
defined status ISTHKK=22; 

• sea qiiarks from interacting target nucleons defined in XKSAMP (number IXTS); 
PHKK(3,...)=PHKK(4,...) contains the actual momentum fraction, VHKK the position of the 'mother' hadron; 
defined status ISTHKK=32; 

• characteristics of the individual parton-parton chains (before hadronization) from subroutines KKEVW, 

KKEVSV, KKEVVS and KKEVSS; for each chain there are three entries: 

(1) two entries for the quark systems forming the chain; 

PHKK gives their 4-momenta; the status of the corresponding quark system is increased by 100 as compared to 
the previous entry from the subroutine XKSAMP (i.e. ISTHKK=121, 122,131 or 132 ,resp.); 

(2) one entry for the complete chain; 

PHKK gives the total 4-momentum, the 'particle' type for chains is defined to be IDHKK=88888; 

the actual status ISTHKK points to the chain generating subroutine: 

ISTHKK=3 for chains from KKEVW (constructed from valence quark systems), 

ISTHKK=4 for chains from subroutine KKEVSV (sea- valence chains), 

ISTHKK=5 for chains from subroutine KKEVVS (valence-sea chains), 

ISTHKK=6 for chains from subroutine KKEVSS (sea-sea chains); 

• hadrons from the hadronization of chains, entries from subroutines HADRSS, HADRVS, HADRSV, HADRW; 
assignment of values to all arrays of /HKKEVT/, status ISTHKK=1; 

• hadrons from resonance decay in JETSET or in subroutine DECHKK 

(presently called after completion of the primary interaction of the projectile treated in KKEVT); 
the status of decaying hadrons is changed to ISTHKK=2, added decay products have ISTHKK=1; 

• hadrons from intranuclear cascade interactions (monitored by FOZOKL): 

the status of interacting secondary is changed to ISTHKK=2; interacting target nucleons get ISTHKK=18, final 
state hadrons have the status ISTHKK=1. 
Particular cases: 
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(i) If a given secondary interaction is found to be forbidden because of the Pauli principle the initial state 
particles are stored in /HKKEVT/ with their original properties, but the actual position; so they may 
participate in further intranuclear interactions. 

(ii) One (or two) nuclcons from a secondary interaction cannot escape from the nuclear potential, but the 
particular collision is not forbidden by Pauli's principle (i.e. several nucleons knocked out of the nucleus 
already): 

Store the nucleon(s) with the actually generated momentum at the collision site, assigning the status 
ISTHKK=15 (16) for interactions in the target (projectile) nucleus. Those nucleons are available as target 
(projectile) nucleons in subsequent steps of the intranuclear cascade development. 

(iii) Negative particles with energies too low to escape from the potential are forced to be absorbed within 
the nucleus (comp. Section 2 and Appendix); absorbed Tr~,K~ and p are characterized by the status 
ISTHKK=19. 

• Evaporation nucleons, nuclear fragments and residual target nuclei from the interface to the evaporation mudule. 

(i) The excited residual nuclei before the evaporation step get ISTHKK(i)=1000 and IDHKK(z)=80000 , the 
mass number A and charge Z of the nucleus are given by IDRES(i) and IDXRES(i). 

(ii) The evaporation protons and neutrons are stored with ISTHKK(i)=-l. 

(iii) The nuclear fragments from evaporation are stored with ISTHKK(i)=-l and IDHKK(j)=80000, the mass 
number A and charge Z of the nucleus are given by IDRES(i) and IDXRES(z). 

(iv) The deexcitation photons are stored with ISTHKK(i)=-l. 

(v) The stable residual nuclei after the evaporation step are stored with ISTHKK(i)=1001 and 
IDHKK(i)=80000, the mass number A and charge Z of the nucleus are given by IDRES(i) and IDXRES(i). 

The information from this common block allows a rather dctaik^d reconstruction of the actual event history, and is 
particularly useful for consistency tests and debugging. An example is discussed in Subsection B.3. 
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B.3: Sample event history from the common block /HKKEVT/ 

In this subsection we discuss a typical event history as stored in the common block /HKKEVT/ according to the 
conventions discussed in the preceding subsection B.2. The corresponding event was sampled for the interaction of 
a 200 GeV proton with a Sulfur nucleus. The following quantities from the common are listed for each entry; the 
number i of the entry, its status ISTHKK(i) and type IDHKK(i), the pointers JM0HKK(2,i) and JDAHKK(2,i) 
to the 'mother' and 'daughter' particles, resp., and the array PHKK(5,i), for final state particles containing the 4- 
momentum and their mass. The last 4 integers are IDRES(ij, IDXRES(i), NOBAM(z) and IDBAM(i). 
The first 33 entries are made in the subroutine KKEVT for the incident proton and the 40 nucleons from the Sulfur 
nucleus. According to our conventions interacting hadrons get the status 11 (projectile) and 12 (target) - there are 
two target nucleons participating in the primary interaction (entries 2 and 25, i.e two neutrons). There are some 
target nucleons which got the status ISTHKK(i)=18, this are target nucleons participating in the formation zone 
intranuclear cascade. 
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The next entries come from the subroutines XKSAMP and FLKSAM (in module DPMNUC3) assigning x values and 
flavor to the quarks participating in the interaction: The projectile proton is split into a d valence quark, an valence 
diquark {uu) (both with status ISTHKK=21), and two quarks d and Jfrom a colorless sea pair (status ISTHKK=31). 

The interacting target nucleons are split in valence quark-diquark systems: the neutron (entry 2) into d-{ud), the 
neutron (entry 25) into d-{ud). The x values given for each quark system as PHKK(3,i)=PHKK(4,i) add up to one 
for each individual hadron. 
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The following entries come from the routines constructing the colorless parton-parton chains. First the chains are 
considered (in subroutine KKEVSV) which are composed from projectile sea quarks and target valence quarks. In 
this example there are two such chains: d-{ud) and d-u. The entries are made subsequently for the partons forming 
the chain and the constructed chain itself, with the array PHKK(5,i) containing the corresponding 4-momenta and 
the mass, resp. The 4-momenta are defined in the cms system of the incident hadron and a single nucleon (without 
Fermi momentum). By definition these chains get the status ISTHKK=4, the status of quarks/quark systems at chain 
ends is increased by 100 compared to their previous entries. 

Please note: quarks, diquarks and chains have no BAM JET code, therefore the last entry in the line (IDBAM(i)) 
picks up values without significance 
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In the same manner the next entries describe the chains constructed from valence quark systems of both the projectile 

and a target nucleon (neutron in this example) in the subroutine KKEVVV. Entries for all chains are marked by the 
identity IDHKK=88888, valence-valence chains get the status ISTHKK=3. 
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The next entries arc generated by the routines HADRxx monitoring the chain decay. 

First the decay of the (two) sea-valence chains is handled via HADRSV. With ISTHKK(i) = l we find the final stable 
hadrons, with ISTHKK(i)=2 we find the resonances, which have decayed already within the JETSET fragmentation 
code. The stable hadrons resulting from resonance dacay point to the decaying resonance (IMOHKK(l,i)). Also stable 
hadrons might have ISTHKK(z)=2, in this case they have interacted subsequently in the formation zone cascade. 

After event completion the 4-momenta are given in the laboratory system. 
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Next we find the entries from the formation zone cascade in subroutine FOZOKL. The status of the initial state target 
nucleon is changed to ISTHKK=18. 

In some of the entries we find the status code ISTHKK(i)=16. This are nucleons from secondary interactions not 
able to escape from the nucleus. These nucleons are available as targets for subsequent collisions. 
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The final entries in the COMMON block refer to the evaporation step of the excited residual target nucleus. With 
status code ISTHKK(i)=1000 we find the excited residual nucleus. Evaporation nucleons, deexcitation photons and 
nuclear fragments follow with status code ISTHKK(i)=-l and the stable residual nuclei follows with status code 
ISTHKK(i)=1001. 
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Eventually the actual final state hadrons, photons and nuclei are obtained from the common block /HKKEVT/ by 
picking up all particles with the status ISTHKK=1, ISTHKK=-1 and ISTHKK=1001. 
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APPENDIX C : 
FURTHER OPTIONS FOR INTERNAL USE 



Code word = 'PAULI' 

This option monitors the inclusion of PauU's principle for the intranuclear cascade. 
WHAT(l) : Pauh principle active for 1.0 (default) 

WHAT(2) : triggers test prints for debugging if greater than 0.0; 
default : 0.0, i.e. no additional printout. 

Code word ^ 'INTPT' 

This option monitors the additional sampling of transverse momenta for partons; 
WHAT(l) : 1.0 activates the sampling of parton values (default). 

Code word = 'TECALBAM' 

This option triggers tests of the hadronization routines; additional input cards are required (see 
subroutine TEC ALB), to be given immediately after the code word definition; WHAT parameters 
have no meaning. 

Code word = 'HADRIN' 

Via this option the secondary interactions generated from FHAD may be forced to be all inelas- 
tic/elastic ones. 

WHAT(l) = : inelastic/elastic as monitored by FOZOKL 

1 : inelastic collisions forced 

2 : elastic collisions forced. 

Code word = 'OUTLEVEL' 

This option monitors the printout of intermediate results for diagnostics; 

WHAT(I) > 1: diagnostics for remaining minor problems is printed (mainly consistency/ accuracy 
problems in kinematical calculations). 

Code word = 'INFOREJE' 

This option monitors the printout of rejection information 
WHAT(I) = 1: rejection information is printed. 

Code word = 'RECOMBIN' 

Combinatoric transformation of s-s and v-v chains into s~v and v-s chains 

WHAT(l) = no recombination, 

What(l) — 1 recombination implemented. 
Default: 0. 

Code word = 'CMHISTO' 

The events in COMMON HKKEVT are sampled in the CMS frame 

WHAT(l) = sampling in lab frame, 
What(l) = 1 sampling in CMS frame 
Default: 0. 

Note: If WHAT(1)=1 there is no possibility to treat the formation zone cascade. See cards TAUFOR 
and PROJKASK 

Code word = 'FLUCTUAT' 

implement cross section fluctuation (see Ref. p3[) 

WHAT(l) = no fluctuations, 

What(l) = 1 fluctuations implemented. 
Default: 0. 
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APPENDIX D : 

FURTHER OPTIONS FROM THE DTUJET-99 CODE WITHIN DPMJET 



Code word: STRUCFUN 

Defines the structure functions used in tlie sampling of hard constituent scattering ||2j 

WHAT(l): = ISTRUF 

ISTRUF = 21 : Glueck,Reya,Vogt with K = 1: GRV94LO 

ISTRUF = 22 : Glueck,Reya,Vogt with K = 2: GRV98LO 

ISTRUF 23 : CTEQ Collaboration (1996): CTEQ96 

ISTRUF = 221,222,223 : as above with energy dependent p± threshold value 

The default is 222 . Only for this structure function we have a new fit. Therefore the use of 

STRUCFUN with ISTRUF diflerent from 222 is not recommended in DPMJET-II.5. 

An energy dependent cutoff avoids a hard scattering cross section too large to be treated in our 
simple eikonal approximation. To use this new option the number 200 has to be added to the chosen 
ISTRUF value. 

Code word: PSHOWER 



This card determines whether hard partons initiate showers in connection with JETSET fragmenta- 
tion. 

WHAT(l): = IPSHOW 

IPSHOW : Generation of hard parton showers suppressed. 
IPSHOW = 1 :Hard parton showers are included. 
The default is 1. 

Code word: GLUSPLIT 

Prevents splitting of gluons into quark-antiquark pairs. 

WHAT(l): = NUGLUU, Default: 1 

NUGLUU=1 :Only one jet in hard gluon scattering. 

Recommended: NUGLUU=1. 
WHAT(2): = NSGLUU, Default: 

NSGLUU=0 :Two jets in soft sea gluon jets. 

NSGLUU=1 :Only one jet in soft sea gluon jets. 

Recommended: NSGLUU^O. 

Code word: SIGMAPOM 

Defines options and/or demands a test run for the calculation of the DTU model unitarizing the soft 
and hard hadronic cross sections. The test run is for calculating the total and inelastic cross sections 
as functions of the collision energy as well as initializing and testing the sampling of multi-Pomeron 
events at some typical energies. Without the test run only the initialization at the energy defined in 
the run is done. 

WHAT(l): = ITEST, test run for ITEST = 1. 

WHAT(2) = ISIG, characterizes the soft and hard input cross sections for the unitarization (see 
SUBROUTINE SIGSHD). 

Recommended and default and only remaining option in DPMJET-II.5 : ISIG = 10 
WHAT(3) = IPIM characterizes the method to calculate the distribution in the number of soft 
and hard ri/j Pomerons 

Default and recommended and only remaining option in DPMJET-II.5: IPIM = 482 
WHAT (4) = LMAX Default: 30 

LMAX :Maximum considered number of soft Pomerons (Less than 26). 
WHAT(5) = MMAX Default: 100 

MMAX :Maximum considered number of hard Pomerons (Less than 101). 
WHAT(6) = NMAX Default: 2 

NMAX :Maximum considered number of triple-Pomerons (Less than 13). 
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